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ABSTRACT 

The relationships between solar flare parameters (total importance, time duration, flare index, 
and flux) and sunspot activity (R z ) as well as those between geomagnetic activity (aa index) and 
the flare parameters can be well described by an integral response model with the response time 
scales of about eight and thirteen months, respectively. Compared with linear relationships, the 
correlation coefficients of the flare parameters with i? z , of aa with the flare parameters, and of aa 
with R z based on this model have increased about 6%, 17%, and 47% on average, respectively. 
The time delays between the flare parameters with respect to R z , aa to the flare parameters, 
and aa to R z at their peaks in solar cycle can be predicted in part by this model (82%, 47%, 
and 78%, respectively). These results may be further improved when using a cosine filter with 
a wider window. It implies that solar flares are related to the accumulation of solar magnetic 
energy in the past through a time decay factor. The above results may help to understand the 
mechanism of solar flares and to improve the solar flare prediction. 

Subject headings: solar physics; solar activity; sun spots; flares; geomagnetic activity 



Introduction 



Solar flares are powe r ful e ru ptions of solar 



activ ity 
1994 



( Ozgiic fc Atacl Il989l ; (Mikic fe Linker 



Jain et all |2010|; iFand 120111 ) occurring 
on time scales of minutes up to a few hours 



([Chandra et al.ll201l[ ) and may produce a series of 
solar-terrestrial effects, which may be hazardous 
to both spacecraft and astronauts. Understand- 
ing the mechanism of solar flares and forecasting 
them are important for both solar physics and 
geophysics. Several mechanisms have been pro- 
posed to explain the eruptions of solar flares, such 
as the photospheric co nverging and shear motions 
dMikic fc Linkerlll994l). 1 lux emergence and cancel- 
lation (|Gan et ah 1993t Zhang et al.ll2001l) . catas- 
trophe model (|Forbedfl990ir and Kink instabil- 
ity of coronal flux ropes (|Sakurailll976l : iLi fe Gan 
20111 ). The magnetic reconnection plays an impor- 



Fandl2010h . 

To qu antify the daily flare activity over 24 hours 
per day, iKleczekl (|l952t) introduced the 'flare in- 
dex' defined as 

Q = ixt, (1) 

where represents the intensity scale of impor- 
tance and the duration (i n minutes) of the flar e 
(jKnoska fc Petrasekl 11984 lAtac fc Ozgficj 1 1998ft ■ 



This relationship is assumed to g ive roughly th e 



total energy emitted by a flare (|Kleczekl [l952ft . 



The solar flare activity is found to be c losely 
correlated with sunspot s (jOzgiic fc Atad Il989t 



Feminella fc Storinil 119971 ). Larger flares appear 
often near larger and more complex active re- 



gions (lMcIntoshlll990l; iBachmann fc White! fl994 ; 



tant role in triggering solar flares (ILin et al. 
Wheatland fc Litvinenko 2001: Forbes et al 



1995 



2006; 



Norauistl l201ll) . Sunspot activity is a striking 
manifestation of magnetic fields on the Sun, as- 
sociated with the main sites of solar-activity phe- 
nomena (|Moradi et all |2010| ) and related to the 



energy supplied into the corona (|de Toma et al 
2000; lTemmer et al . 1 120031 ) . Studying the relation 



1 



ship between solar flares and sunspot activity is 
useful to understand and predict the former. The 
flare frequency of occurrence is often predicted 
by sunspot groups or numbers ( Mcintosh 1990; 



2002; Cui et al. 2006; Yu et al 



Huang et al.l 120101 )" for increasing applica- 



Galla gher et al 
2009; 

tions in space weather 

Solar activity is well know n to be at the ori- 



gin of geomagnetic activity (jSnvder et al.l [1963; 



Crooker et al.l 119771 ). Studying the relationship 



between solar activity, as represented by the Inter- 
national sunspot number (R z ), and geo magnetic 
activity, as represented by the aa index ( Mayaudl 
19721 ). is useful for understanding the formation 
of t he latter and the mechanism of solar cy- 
cle (iFevnman fc Crooker 1978; iLegrand fc Simonl 



1989UDvJl2011aHDu fc Wandl2010U2011aM . Con 

ventionally, the relationship between aa and R z 
is often analyzed by point-point correspondence. 
However, some questions are hardly understood 
such as the significant incr ease in the aa index 
over the twentieth centur y (Fcynman & Crooker 
19781 IClilverd et alJll998t iLukianova et al.ll200gF 7 
and the yariations in the correlation between aa 
and R z (Bor eho-Filisetti et alJll992t Echer et al 



2004t lDull2011bl ). It is found that these phenom- 



ena can be well explained by an in t egral r esponse 
model recently presented by IDuI ( 2011c ). The 
value of aa depends not only on the present R z 
but also on past values. 

The geomagnetic activity results from various 
phenomena which are rela ted to the interplane - 



tary magneti c field (IMF. IStamper et al.l 11 999) 



solar wind (Svalgaardl 19771 ; ILegrand fe Si 



19891 : iTsurutani et al.lll995l). Co ronal Mass Ej 



tion (CME. ILegrand fc Simonl Il989l) . galactic 
cosmic rays ( Stamper et ah fll999h. and others 



cosmic rays (Stamper et ai.i liiJHiJl). ana otner s 
( Leerand fc Simon! 1989; Stamper et al.l Il999h . 
Gosling! (|l993l ) pointed out that CMEs, rather 
than flares, were the critical element for large ge- 
omagnetic storms, interplanetary shocks, and ma- 
jor solar ener getic particle (SEP) event s, which 
was arg ued bv (Richardson fc Canel (|2002l) . 

This study analyzes the relationships between 
solar flare parameters (Section [2]) and R z as well 
as the relationships between the aa index and the 
flare param eters using an integral response model 
( Du 2011cl ) in Sections I3.HI3.4I Conclusions are 
summarized finally in Section 01 
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Fig. 1— (a) I (solid) and T (dotted) since March 
1976, with a correlation coefficient of r — 0.89. (b) 
Q (solid) since July 1966 and F (dotted) since July 
1997, with a correlation coefficient of r = 0.67. (c) 
R z (solid) and aa (dotted) since July 1966, with a 
correlation coefficient of r — 0.40. 

2. Data 

The data used are the tim e series of mo nthly 
mean geomagnetic aa indesQ (|Mavaudlll972l) . the 
international sunspot number (R z and solar 
flare parameters based on Geostationary Opera- 
tional Environmental Satellite (GOES) soft X-ray 
flares shown as follow^. 

(i) /: total importance of flares, / = 100X + 
10M + C + 0. IB where X, M, C and B arc 
the flare classes ( Cui et aUl2006l) . 



(ii) T: time duration of flare (in minutes) . 

(iii) Q: the ' flare index' from Equation ([1]) by 



Kleczekl (|1952| ) 



(iv) F: flux from event start to end (in J/m 2 ) . 

These parameters are first summed over each 
day and, then, averaged over each month to ob- 
tain the monthly means of the daily integrated 
quantities. To filter out high frequency variations 
in the data, the parameters are smoothed with 
the commonly used 13- month running mean tech- 
nique. The solar flare parameters since July 1966 



1 ftp://ftp. ngdc.noaa.gov/STP/SOLAR_D ATA/RELATED- 
JNDICES/AAJNDEX/ 
2 http: / / www. ngdc . noaa. gov /stp/ spaceweather . html 
3 ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/ 
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are shown in Fig.QJa) and (b). It is seen that these 
parameters are well correlated. For example, / is 
well correlated with T (r = 0.89, Fig. Ufa)), and 
Q is well correlated with F (r = 0.67, Fig.ETJb)), 
both being significant at the 99% level of confi- 
dence. Figure QJc) depicts the time series of R z 
(solid) and aa (dotted) with a correlation coeffi- 
cient of r — 0.40. 



3. Results 

It is well known that solar flares tend to 
l ag behind sunspot activ i ty by s everal months 
(IWheatland fc Litvinenko 2001; Temmer et al 



2003 ) or even a few years ( Wagnerll988t Aschwanden 
1994T ). To have a better understanding of the re- 
lationships and time delays between solar flares 
and R z , we employ the following integral response 
model joi] 1201 let) to study the relationships be- 
tween the flare parameters (P = I, T, Q and F) 
and i? z , 



y(t) = Djt^xWe-WTtf + yo 



(2) 



where yo is a constant, reflecting the part of y = P 
that is uncorrelated to x = R z (related to other 
phenomena); D is the 'dynamic response factor' of 
y to x, representing the initial generation efficiency 
of y by x (dy /dx\ t >=t)] and r is the 'response time 
scale' of y to x, indicating the dependence of the 
current y(t) on the past x(t') through a time de- 
cay factor e~(* - * ^ T (r = reflects the point-point 
correspondence of y to x, i.e., the current y(t) is 
only related to the current x(t); r = +oo repre- 
sents that y is uncorrelated to x). In application, 
both y and x are discrete variables. Therefore, 
we use the second formula in Equation @ with 
the summation being taken over from the starting 
time (to) of the series (see Fig. [J) to time t. The 
three parameters (D, r and yo) are determined 
by a nonlinear least-square fitting algorithm. Be- 
sides, as the geomagnetic activity (aa index) often 
lags behind solar flares by several months, the re- 
lationships between aa and the flare parameters 
are also analyzed by the same model. 

3.1. Relationship between R z -I-aa 

First, we analyze the relationship between y = 
I and x — R z since March 1976 (to) with Equation 
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Fig. 2. — (a) I (solid), i? z (dashed), and the re- 
constructed series (If, dotted) by Equation ((3]). 
The correlation coefficients of I with R z and If 
are r = 0.87 and n = 0.88, respectively. The lag 
times of I (If) to R z at their peaks for Cycles 21-23 
are L x = 31, 24, 17 (L n = 4, 8, 5) months, (b) Sim- 
ilar results for the relationship between aa (solid) 
and I (dashed). The correlation coefficients of aa 
with / and the reconstructed series aaf by Equa- 
tion (U]) are r = 0.61 and rf — 0.74, respectively. 
The lag times of aa (aaf) to I at their peaks for Cy- 
cles 21-23 are L 2 = 4, 2, 23 (L i2 = 5, 4, 8) months, 
(c) For the relationship between aa (solid) and 
R z (dashed). The correlation coefficients of aa 
with R z and the reconstructed series aaf by Equa- 
tion ([5]) are r = 0.51 and rf = 0.73, respectively. 
The lag times of aa (aaf) to R z at their peaks for 
Cycles 21-23 are L = 35,26,40 (L { = 27,29,28) 
months. 

@ in the form of 



I(t) 



AEL^(i')^ (i - t ' )/ HI„. (3) 



Figure [2{a) plots the reconstructed series If (dot- 
ted) of I (solid) from R z (dashed) by Equation ((3|) . 
Although the correlation coefficient between I and 
Rz (fQ = 0.87) has not been significantly improved 
by this model (rf = 0.88), the lag times of I to 
R z at their peaks (time differences between the 
peak timings )_for Cycles 21-23 (Li = 31,24,17 
with a mean L\ = 24 months) can be predicted 
in part by Equation (J3|) as shown in Fig. [U[a) for 
the corresponding ones in brackets (Lfi = 4,8,5 
with a mean Lf\ =6). It implies that the cur- 
rent flares are related to the accumulation of so- 
lar magnetic energy in the past through a time 
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decay factor. Active magnetic structures may 
evolve from the photosphere to uppe r chromo 



sphere with different speeds and tim es (|Lin et al 
19951; IWheatland fc Litvinenkoll200l"h . 



The relationship between y 
be fitted by 



aa and x — I can 



aa(t) 



as shown in Fig. [UJb): aa (solid), / (dashed), 
and the reconstructed series aaf (dotted) by Equa- 
tion (j4j. One can see that aaf reflects well the 
profile of aa. The correlation coefficient between 
aa and aaf (rf = 0.74) is higher than that be- 
tween aa and / (tq = 0.61). About half of the lag 
times of aa to / at their peaks for Cycles 21-23 
(L2 = 4, 2, 23 with a mean L2 = 10) can be pre- 
dicted by Equation (j4|) as shown in Fig. EJb) for 
the corresponding ones in brackets (Lf2 — 5,4,8 
with a mean Lf2 = 6). 

The relationship between y = aa and x = R z is 
analyzed by using the following equation, 



»(*) = DJ2 t t , =t0 R z (t')e-^ t '^ + 



aao- 



(5) 



Figure [^c) illustrates aa (solid), R z (dashed), and 
the reconstructed series aaf (dotted) by this equa- 
tion. The correlation coefficient between aa and 
aaf (rf = 0.73) is much higher than that between 
aa and R z (r — 0.51). The lag times of aa to R z 
at their peaks for Cycles 21-23 (L = 35, 26, 40 with 
a mean L = 34) can be well predicted by Equa- 
tion ([5]) as shown in Fig. [Jfc) for the correspond- 
ing ones in brackets (Lf = 27, 29, 28 with a mean 
Li = 28). The above results are listed in Table Q] 
in which a refers to the standard deviation, the 
last column indicates the relevant averages of fit- 
ted/observed lag times at the corresponding peaks 
over Cycles 21-23 (L[/L), and the last three rows 
represent the relevant averages of the parameters 
for the relationships between P-R z , aa-P, and aa- 
R Zl respectively, where P = I,T, Q, F. 

3.2. Relationship between R z -T-aa 

The relationship between R z -T-aa since March 
1976 (t ) can also be analyzed by the technique 
in the previous section, with the results shown in 
Fig. |3] and Table [TJ The following can be noted. 

(i) The correlation coefficient of T with the re- 
constructed series Tf (rf = 0.89) from R z by 
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Fig. 3. — Similar to Fig. [2] for the relationship 
between R z -T-aa 
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Fig. 4. — Similar to Fig. [5] for the relationship 
between R z -Q-aa. 

Equation @ is higher than that of T with 
R z (r a = 0.79). 

(ii) The correlation coefficient of aa with the re- 
constructed series aaf (rf = 0.74) from T by 
Equation @ is higher than that of aa with 
T(r = 0.66). 

(iii) The lag times of T to R z (Lfi/Li = 25/21), 
aa to T (L12/L2 = 5/12), and aa to i? z 
(Lf/L — 28/34) at their peaks for Cycles 
21-23 can be approximately predicted by the 
model. 

3.3. Relationship between R z -Q-aa 

Figure 0] illustrates the relationship between 
R z -Q-aa since July 1966 (to) by using the tech- 
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Table 1: Fitted Results of the Integral Response Model for the Flare Parameter: P = J,T, Q, F. 
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b Average of the corresponding parameters for the relationships between P (= I,T, Q,F) and R z . 



niquc in Section [3. II One can note the following. 



(i) 



The correlation coefficient of Q with the re- 
constructed series Qf (r f = 0.90) from R z 
by Equation ([2]) has not improved in com- 
parison to that of Q with i? z (tq = 0.90), 
implying t hat Q and R z peak nearly at the 
same time ( Kleczeklll952l ). 



(ii) The correlation coefficient of aa with the re- 
constructed series aa{ (rf = 0.58) from Q by 
Equation ([2]) is much higher than that of aa 
with Q (r = 0.37). 

(iii) The correlation coefficient of aa with the re- 
constructed series aaf (rf = 0.64) from R z 
by Equation ([2]) is much higher than that of 
aa with R z (tq = 0.32). 

(iv) The lag times of Q to R z (i fl /ii = 0/5), 
aa to Q (L{ 2 /L 2 = 24/29), and aa to R z 
(Lf/L = 31/34) at their peaks for Cycles 
21-23 can be predicted in part by Equation 
©• 

3.4. Relationship between R z -F-aa 

Figure [5] shows the relationship between Rz-F- 
aa since July 1997 (to) using the technique in Sec- 
tion [3J] One sees the following. 



(i) The correlation coefficient of F with the re- 
constructed series F[ (rf = 0.89) from R z by 
Equation ^ is higher than that of F with 
R z (r = 0.79). 

(ii) The correlation coefficient of aa with the re- 
constructed series aaf (rf = 0.75) from F by 
Equation ([2]) is equal to that of aa with F 
(ro = 0-75). 

(iii) The correlation coefficient of aa with the re- 
constructed series aaf (rf = 0.79) from R z 
by Equation @ is much higher than that of 
aa with R z (r = 0.62). 

(iv) The lag times of F to R z (ifi/ii = 25/17), 
and aa to R z (Lf/L = 26/40) at their peaks 
for Cycle 23 can be predicted in part by 
Equation @. While the lag time of aa to 
F (L n /L 2 = 0/23) at their peaks for Cycle 
23 has not been predicted by Equation (J2]) 
due to the great fluctuations in both aa and 
F. 

These results imply that solar flares depend 
not only on the present but also on past solar 
activities (i? z ), reflecting the long-term evolution 
characteristics of solar magnetic field structures 
(energy) evolvi ng from the photosphere to upper 
chromosphere (|Donnellvl 119871 : IZhang et al.l 12007 ; 
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Fig. 5. — Similar to Fig. [2] for the relationship 
between R z -F-aa. 
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Fig. 6. — Same to Fig. [2] but using the cosine filter 
(Equation ©). 



Lin|[2009). The correlations between solar flare pa- 
ramete rs and R z are not simply d ue to the time 
shifts (jBachmann fc White] Il994l ). Solar flares 
may play a role for the formation of geomagnetic 
activity (aa) from solar (magnetic field) activity 
(R z ), although the processes are not completely 
clear (jCliver fc Hudsonll2002l ). 



4. Discussions and Conclusions 

In this study, we investigated the relationships 
between the solar flare parameters (P = I, T, Q 
and F) and sunspot activity (i? z ), and between ge- 
omagnetic activity (aa) and the flare parameters 
via the integral response model (Equation ([2])). 
The results indicate that (i) the correlation coef- 
ficients between the flare parameters and R z have 
increased about 6% from To = 0.84 to Tf = 0.89 
on average when using Equation @ and the time 
delays at their peaks for Cycles 21-23 can be well 
predicted by this model, Ln/Li = 14/17 = 82%; 
(ii) the correlation coefficients between aa and 
the flare parameters have increased about 17% 
from To — 0.60 to Tf — 0.70 on average when 
using Equation ([2]) and half of the time delays 
at their peaks can be predicted by this model, 
Lf 2 /L 2 = 9/19 = 47%; and (iii) the correlation 
coefficient between aa and R z has increased about 
47% from To = 0.49 to Tf — 0.72 on average 
when using Equation @ and the time delays at 
their peaks can be well predicted by this model, 
L { /L = 28/36 = 78%. This model might be used 
to improve the solar flare prediction, which should 



be studied in future. 

It is seen in Fig.[5Ja) and TableQ]that the time 
delays between / and R z at their peaks for Cycles 
21-23 have not been well predicted by the model 
(6/24 = 25%). This is due to the large fluctuations 
in the data. To suppress further the fluctuations, 
we introduce a cosine filter with the weights given 

by 

W c (At) = ^cosA (6) 



■ib 



2!> 



for b = 24 months. Since f_ b ^ cos(^) = 1, the 
weights Wc(At) are normalized. Using the series 
smoothed by this filter, we re- analyze the results 
in Fig. [21 as shown in Fig. [Bl 

The time delays between I and R z at their 
peaks are now better predicted, Ln/Li = 5/15, 
5/1 and 4/3 for Cycles _21-2^ (Fig. E^a)), respec- 
tively, with a mean of L lf /L 1 = 4.7/6.3 = 74% 
which is much higher than the original one (6/24 = 
25%). The time delays between aa and / at their 
peaks are also better predicted, L21/L2 = 11/16, 
13/13 and 19/27 for Cycles_21-23 (Fig.^b)), re- 
spectively, with a mean of L^ijL^ = 14.3/18.7 = 
77% which is higher than the original one (6/10 = 
60%). In Fig. [6jc) , the time delays between aa and 
R z at their peaks are predicted as L[/L = 19/31, 
17/14 and 16/30_for Cycles 21-23, respectively, 
with a mean of Lf/L = 17.3/25 = 69% which is 
smaller than the original one (28/34 = 82%) due 
to the great lag time of aa to R z (about 30 months) 
and other sources of aa. 

In Equation @, the output y depends on the 
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past values of input x (r > 0) rather than only 
the current value (r = 0). The stronger the 
input (x), the more it contributes to the out- 
put (y), an d the longer the lag time of y to x 
(flhH 1-201 id). Therefore, solar flares are related 
to the accumulation of solar magnetic energy in 
the past ra ther than the simple time shifts of oc- 
currences ( Bachmann fc White 1994). The av- 
erage response time scale of flare parameters to 
R z in this model (r = 8) is close to the coro- 
nal response time (~ 10 months) derived from a 
model for dynam ical energy balance in the flar- 
ing solar corona (IWheatland fe Litvinenkd 12001 



Litvinenko fc Wheatlandll2004l ) 

There are various types of active regions in a 
solar cycle. Small active regions of simple mag- 
netic structure are short-lived and produce mi- 
nor solar flares, while large active regions of com- 
plex magnetic structure are long-lived and pro- 
duce major solar flares (and hard X-ray flares). It 
is seen in Fig. [lja) that T is well correlated with / 
(r = 0.89), with the regression equation given by 



T = 67.9 ± 2.2 + (1.94 ± 0.05)/. 



(7) 



According to the above discussions, minor (low- 
energy) solar flares lag behind the input R z shorter 
times with shorter durations while major solar 
flares lag behind R z longer times with longer du- 
rations. Therefore, (i) the time delays between 
flare activities and sunspot activity come mainly 
from the major flares rather than the weak ones; 
(ii) major flares tend to have longer durations 
and may occur u ntil quite late in the decay phas e 
of a solar cycle (|Temmer et all 120031 : lTan1l201l[ k 



and ( hi) the upper chromospheric activity in- 
dices (|Donnellvlll987t iBachmann fe White! Il994l) 
and the solar flares (jWheatland fc Litvinenko] 



20011 : iTemmer et all 120031 ) tend to lag behind the 



sunspot number by several months in a hierarchy 
manner (|Bachmann fc White! Il99l . 



Although it is unc lear how solar flares affect ge- 
omag netic activities ( Goslindl993 ; Oliver fc Hudson 
l2002h . it is apparent that geomagnetic activities 
are well correlated with the solar flares. For ex- 
ample, aa is well correlated with F (ro = 0.75). 
As flares are unable to travel to 1 A.U., streams 
of matter emanating from large flares were con- 
sidered as the prime cause of geomagnetic storm s 
(lHaldll93lMChapmadll950l:lPudovkin et al.lll977l) . 
However. iGoslind (|l993l) argued that CMEs, not 



flares, were the critical element for large geomag- 
netic storms, interplanetary shocks, and major so- 
lar energetic particle (SEP) events. In fact, solar 
flares may affect geomagnetic activities via dif- 
ferent processes related to the flare brightening, 
eruptin g, particle ejections, an d other unknown 
effects (|Cliver fc Hudson! 120021 ). Therefore, the 
relationships between geomagnetic activity (aa) 
and solar flares can also be well described by 
Equation ©. Since geomagnetic activity (aa) 
can be resulted from var i ous activity phenomen a 
(jLegrand fc Simon! Il989l : iTsurutani etail Il995l) . 
the geomagnetic activity is the integral of the 
effects of all these phenomena, including solar 
winds, CMEs, solar flares and others. The lag 
time of aa to solar flare has not been well pre- 
dicted by the model (9/19) due to the additional 
effects of other activities. While the lag times of 
both solar flare and aa to R z at their peaks have 
been well predicted by the model (14/17, 28/36) 
because the solar magnetic field activity is the 
main source of them. 

The main conclusions can be drawn as follows, 

(i) The relationships between the flare param- 
eters (P — I,T,Q,F) and sunspot activ- 
ity (R z ) can be well described by an inte- 
gral response model (r = 0.89) with a mean 
response time scale of about eight months. 
The time delays between the flare parame- 
ters and R z at their peaks can be well pre- 
dicted by this model (82%). 

(ii) The relationships between geomagnetic ac- 
tivity (aa) and the flare parameters can be 
better described by this model (r = 0.70) 
with a mean response time scale of thir- 
teen months than by a linear dependence 
(r = 0.60). The time delays between aa and 
the flare parameters at their peaks can be 
predicted in half by this model (47%). 

(iii) The relationship between aa and R z can be 
much better described by this model (r = 
0.72) with a mean response time scale of 
about twenty-six months than by a linear 
dependence (r — 0.49). The time delay be- 
tween aa and R z at their peaks can be pre- 
dicted in part by this model (78%). 



7 



Acknowledgments 

The authors are grateful to the anonymous ref- 
eree for constructive comments. This work is sup- 
ported by National Natural Science Foundation of 
China (NSFC) through grants 10973020, 40890161 
and 10733020, and National Basic Research Pro- 
gram of China (973 Program) through grant No. 
2011CB811406. 

REFERENCES 

Aschwanden, M. J. 1994, Sol. Phys., 152, 53 

Atac, T., & Ozgiig, A. 1998, Sol. Phys., 180, 397 

Bachmann, K. T., & White, O. R. 1994, 
Sol. Phys., 150, 347 

Borello-Filisetti, O., Mussino, V., Parisi, M., & 
Storini, M. 1992, Ann. Geophys., 10, 668 

Chandra, R., Schmieder, B., Mandrini, C. H., et 
al. 2011, Sol. Phys., 269, 83 

Chapman, S. 1950, J. Geophys. Res., 55, 361 

Clilvcrd, E. W., Boriakoff, V., & Feynman, J. 
1998, Geophys. Res. Lett., 25, 1035 

Cliver, E. W., & Hudson, H. S. 2002, J. Atmos. 
Sol. Terr. Phys., 64, 231 

Crooker, N. U., Feynman, J., & Gosling, J. T. 
1977, J. Geophys. Res., 82, 1933 

Cui, Y., Li, R., Zhang, L., et al. 2006, Sol. Phys., 
237, 45 

De Toma, C, White, O. R., & Harvey, K. L. 2000, 
ApJ, 659, 801 

Donnelly, R. F. 1987, Sol. Phys., 109, 37 

Du, Z. L. 2011a, Ann. Geophys., 29, 1331 

Du, Z. L. 2011b, Ann. Geophys., 29, 1341 

Du, Z. L. 2011c, Ann. Geophys., 29, 1005 

Du, Z. L., Wang, H. N. 2010, RAA (Res. Astron. 
Astrophys.), 10, 950 

Du, Z. L., Wang, H. N. 2011a, Sci. China Scr. G- 
Phys Mech Astron, 54, 172 

Du, Z. L., Wang, H. N. 2011b, RAA (Res. Astron. 
Astrophys.), 11, 1482 



Fang, C, Chen, P.-F., Jiang, R.-L., Tang, Y.-H. 
2010, RAA (Res. Astron. Astrophys.), 10, 83 

Fang, C. 2011, RAA (Res. Astron. Astrophys.), 
11, 1377 

Echer, E., Gonzalez, W. D., Gonzalez, A. L. C, 
et al. 2004, J. Atmos. Sol. Terr. Phys., 66, 1019 

Feminella, F., & Storini, M. 1997, A&A, 322, 311 

Feynman, J., & Crooker, N. U. 1978, Nature, 275, 
626 

Forbes, T. G. 1990, J. Geophys. Res., 95, 11919 

Forbes, T. G., Linker, J. A., Chen, J. et al. 2006, 
Space Sci. Rev., 123, 251 

Gallagher, P. T., Moon, Y. J., & Wang, H. 2002, 
Sol. Phys., 209, 171 

Gan, W. Q., Rieger, E., & Fang, C. 1993, ApJ, 
416, 886 

Gosling, J. T. 1993, J. Geophys. Res., 98, 18937 

Hale, G. E. 1931, ApJ, 73, 379 

Huang, X., Yu, D. R., Hu, Q. H, et al. 2010, 
Sol. Phys., 263, 175 

Jain, R., Aggarwal, M., & Kulkarni, P. 2010, RAA 
(Res. Astron. Astrophys.), 10, 473 

Kleczek, J. 1952, Ionospheric Disturbances and 
Flares in the 11-years cycle, Bulletin of the As- 
tronomical Institutes of Czechoslovakia, 3, 52 

Knoska, S., & Petrasek, J. 1984, Contr. Astron. 
Obs. Skalnate Pleso, 12, 165 

Kusano, K., Maeshiro, T., Yokoyama, T., & Saku- 
rai, T. 2004, ApJ, 610, 537 

Legrand, J. P., & Simon, P. A. 1989, Ann. Geo- 
phys., 7, 565 

Li, Y. P., & Gan, W. Q. 2011, Sol. Phys., 269, 59 

Lin, J. 2009, Sci. China, Phys. Mech. Astron., 52, 
1646 

Lin, J., Forbes, T. G., Priest, E. R., & Bungey, T. 
N. 1995, Sol. Phys., 159, 275 

Litvinenko, Y. E., & Wheatland, M. S. 2004, 
Sol. Phys., 219, 265 



8 



Lukianova, R., Alekseev, G., & Mursula, K. 2009, 
J. Geophys. Res., 114, A02105 

Mayaud, P. N. 1972, J. Geophys. Res., 77, 6870 

Mcintosh, P. S. 1990, Sol. Phys., 125, 251 

Mikic, Z., & Linker, J. A. 1994, ApJ, 430, 898 

Moradi, H., Baldner, C., Birch, A. C, et al. 2010, 
Sol. Phys., 267, 1 

Norquist, D. C. 2011, Sol. Phys., 269, 111 

Ozgiic, A., & Atac, T. 1989, Sol. Phys., 123, 357 

Pudovkin, M. I., Zaitscva, S. A., Oleferenko, I. P., 
& Chertkov, A. D. 1977, Sol. Phys., 54, 155 

Richardson, I. G., & Cane, H. V. 2002, J. Geo- 
phys. Res., 107, 1187 

Sakurai, T. 1976, PASJ, 28, 177 

Snyder, C. W., Neugebauer, M., & Rao, U. R. 
1963, J. Geophys. Res., 68, 6361 

Stamper, R., Lockwood, M., Wild, M. N., & 
Clark, T. D. G. 1999, J. Geophys. Res., 104, 
28325 

Svalgaard, L. 1977, in Coronal Holes and High 
Speed Wind Streams, ed. J. B. Zirker, Colorado 
Ass. U. Press, Boulder, P.371 

Tan, B. 2011, Ap&SS, 332, 65 

Temmcr, M., Veronig, A., & Hanslmeicr, A. 2003, 
Sol. Phys., 215, 111 

Tsurutani, B. T. Gonzalez, E. D., Gonzalez, A. L. 
C. et al. 1995, J. Geophys. Res., 21, 717 

Wagner, W. J. 1988, Adv. Space Res., 8, 67 

Wheatland, M. S., & Litvinenko, Y. E. 2001, ApJ, 
557, 332 

Yu, D. R., Huang, X., Wang, H. N., & Cui, Y. M. 
2009, Sol. Phys., 255, 91 

Zhang, J., Wang, J., Deng, Y., & Wu, D. 2001, 
ApJ, 548, L99 

Zhang, L. Y., Wang, H. N., Du, Z. L., et al. 2007, 
A&A, 471, 711 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



9 



